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Power and spectral index anisotropy of the entire inertial 
range of turbulence in the fast solar wind. 
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ABSTRACT 



We measure the power and spectral index anisotropy of high speed solar wind 
turbulence from scales larger than the outer scale down to the ion gyroscale, thus 
covering the entire inertial range. We show that the power and spectral indices at the 
outer scale of turbulence are approximately isotropic. The turbulent cascade causes the 
power anisotropy at smaller scales manifested by anisotropic scalings of the spectrum: 
close to fc -5 / 3 across and k~ 2 along the local magnetic field, consistent with a critically 
balanced Alfvenic turbulence. By using data at different radial distances from the Sun, 
we show that the width of the inertial range does not change with heliocentric distance 
and explain this by calculating the radial dependence of the ratio of the outer scale 
to the ion gyroscale. At the smallest scales of the inertial range, close to the ion 
gyroscale, we find an enhancement of power parallel to the magnetic field direction 
coincident with a decrease in the perpendicular power. This is most likely related to 
energy injection by ion kinetic modes such as the firehose instability and also marks 
the beginning of the dissipation range of solar wind turbulence. 

Key words: turbulence - plasmas - solar wind. 



1 INTRODUCTION 

The solar magnetic field embed ded in the turbulent so- 
lar wind (|Goldstein et alj 1 19951 ) causes several kinds of 
anisotropy. Different physical processes operate parallel 
and perpendicular to this mean field. The fluctuations 
(with respect to the large-scale mean field) of the par- 
allel and perpendicular co mponents of the magnetic field 
scale differently from larg e dNicol et al.ll2009l) to small scales 
|Alexandrova et al.1 [20081 ; ISahraoui et al.1 l2009h . Many dif- 



ferent solar wind quantities, for example; 



( Milano et al.l 2004 


), the magnetic field correlation function 


(Dasso et al. 2005; 


Osman & Horbury 2007; Wcygand ct al. 


2009]), and velocity components (Chapman & HnatH2007h, 



have been shown to vary anisotropically with respect to the 
mean magnetic field direction. 

Recent advances in the study of the solar wind have 
shown that the local magnetic field direction causes spa- 
tial anisotropy in the statistical properties of the fluctuat- 
ing fields along and across the local mean magnetic field. 
This anisotropy can b e measured in correlat ion functions 



dMatthaeus et al.1 


19941), p 


119961; IChen et al. 


l2010bh 



and variation of the spectral 
index depending on the angle to the local mean field 
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(IHorburv et al.ll2008l ; |Podestall2009l ; IChen et al.ll2010bl) ; see 
IChen et al l (|2010al ) for a discussion of the relationship be- 
tween these manifestations of the anisotropy. These re- 
sults have shown that the measured spacecraft power spec- 
tra of the turbulent magnetic field, B, at different an- 
gles, 6b, to the local mean magnetic field are consistent 
with those predicted by t he phenomenological theory of 
iGoldreich fc Sridharl (|l995f ) based on the conjecture of criti- 
cal balance between linear (Alfven) and nonlinear timescales 
in the turbulent cascade. Thus far these studies have focused 
on the inertial range, the small-scale boundary of which is 
defined by the ion gyroradius pi, and have not approached 
the large-scale limit, the outer scale, the point at which the 
pow er spectrum fc" rolls over to a flatter index of a ~ — 1 
(e.g. IHorburv et aflll996l h 



Here w e extend the w a velet-b ased local magnetic field 
analysis of IHorburv et"afl (|2008l ) and iPodestal (|2009l ) to 
larger scales, including the outer scale. We investigate the 
properties of the turbulence, in particular the anisotropy, 
across the entire inertial range and demonstrate for the first 
time the development of the turbulent cascade from its orig- 
inal, nearly isotropic, state at low frequencies all the way 
down to to its final, anisotropic, state near the ion gyroscale. 
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WAVELETS AND THE LOCAL MEAN 
MAGNETIC FIELD 



We use the method devised by IHorburv et all (|2008h and 
explained in more detail by iPodestal (|2009h . The Morlet 
wavelet is used to decompose the time series of magnetic field 
observations into powers at given frequencies and scales. 
Neighbouring wavelet scales differ by a factor of 1.6, making 
the Morlet wavelet close to orthogonal and ensuring a good 
coverage of frequency without oversampling. The wavelet 
scale is used as the length over which to calculate the mean 
magnetic field direction (the local mean field method). The 
power P(f, 9b) calculated using the wavelet is then assigned 
to a bin corresponding to an angle 9 b , the angle between the 
mean field and the sampling direction (the direction of the 
spacecraft flight in the solar wind) an d a Fourier frequenc 



/ associated w ith the wavelet scale |Horburv et alj 200 



Podestal 120091 '). We adopt the Taylor hypothesis (|Tavlorl 
19381 ) that the flow speed is much greater than the sound 
and Alfven speeds and so a time series can be considered to 
be a one-dimensional cut through a time stationary plasma. 
Therefore, our observations are nece ssarily of the reduced 
spectrum (jFredricks fc Coroniti|[l973 ). 

We take periods of continuous Ulysses data from 1995 
while the spacecraft is in fast polar sol ar wind. We use on e 
second resolution magnetic field data jBalogh et al.lll992f ). 
Analysing each period produces a power spectrum ranging 
in spacecraft frequency / between 3.3 x 10 -5 and 2.5 x 10 _1 
Hz. We resolve angles 9b between the observed wavevector 
and the local mean magnetic field in 10° bins between 0° and 
180°. In what follows we only plot the bins 0° < 6 B < 90°: 
angles greater than 90° do not occur frequently and so we 
cannot reliably measure the anti-parallel power spectrum. 

Ulysses observations are made at 1 s cadence, however, 
the important physical scale for kinetic plasma physics in the 
solar wind is the proton gyroscale pi. In order to compare 
different periods directly and to cast our results in phys- 
ically relevant units we convert the spacecraft observation 
frequency into a flow-parallel wavenumber k by dividing by 
the average solar wind speed \ V\ and normalise this by pi\ 



kpi = 



2-Kfpi _ 2TYfV2k B T t m t 
\V\ 



e\V\\B\ 



(1) 



where ks is Boltzmann's constant, Ti is the proton temper- 
ature, rrii is the mass of a proton, e is the charge on an 
electron and |B| is the magnetic field strength. 



3 ANISOTROPY OF THE ENTIRE INERTIAL 
RANGE 

First, we analyse a period of fast polar wind from Ulysses 
data between days 100 and 200 of 1995; during this time 
Ulysses moved from a solar latitude of 28° to 79° and from 
a heliocentric distance of 1.38 AU to 1.93 AU. This long 
interval provides an accurate anisotropic power spectrum at 
the lowest frequencies used here; shorter periods can be used 
if angular resolution is not required at such low frequencies. 
Figure [1] shows the trace of the magnetic field power tensor, 
averaged over periods when the solar wind flow is parallel, 
P\\ (0° < 6 B < 10°) and perpendicular, P ± (80° < 6 B < 
90°) to the local magnetic field calculated using wavelets 
and the average windowed Fourier power calculated between 
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Figure 1. Trace of the wavelet and Fourier power spectra of mag- 
netic field observations from Ulysses for the period between days 
100 and 200 of 1995. Frequencies are converted to wavcnumbers 
using the solar wind velocity and normalised to the ion gyroradius 
pi (Equation [T} . Sec Figure [4] for compensated spectra. 



the same frequencies for the same period. At the smallest 
values of kpi the power is isotropic and all three lines lie 
close together with a spectral index of approximately — 1. At 
kpi ps 3 x 10~ 3 , Pi I begins to diverge from the Fourier power 
and P±. The power anisotropy increases as kpi increases 
from this point; P± and the Fourier power follow each other 
very closely and are a factor of five larger than Pit at the 
largest kpi measured. We stress that the use of wavelets to 
analyse the anisotropy of the magnetic field means that the 
magnetic field is not broken into components parallel and 
perpendicular to the field. Thus the terms Pu and P± do not 
refer to components of the field but to the mean trace power 
in the field when the flow past the spacecraft is parallel or 
perpendicular to the mean field. 

In Figure [2] we show the dependence of spectral index a 
on scale by plotting a determined by a least-squares fitting 
in logarithmic space of a straight line to power measure- 
ments. Each fit uses five wavelet power measurements and 
so covers scales separated by a factor of 10.5. The kpi shown 
are the central values for each five point window, the uncer- 
tainties are calculated as the error from the fitting of the line. 
At large scales {kpi = 1.9 x 10 -4 ), the power spectra have 
an isotropic spectral index a ~ — 1. Above kpi ~ 5 x 10~ 4 , 
the spectrum steepens (diamonds in Figure [2}. As kpi in- 
creases, the spectral indices decrease further and their val- 
ues vary with 9b but only slowly with scale (squares and 
triangle s in Figure [2]). These scale s show a s imila r result to 
those of IHorburv etalj (120081 ) and lPodestal (|2009l ) and con- 
firm that in the inertial range fluctuations have a parallel 
spectral index of —2, while the perpendicular spectral in- 
dex lies between the —5/3 of a crit ically balanced Alfvenic 
cascade (|Goldreich fc Sridhar|[l995l ) and the -3/2 predicted 
by the theory of dynamically aligned Alfvenic turbulence 
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Figure 2. Dependence of the spectral index in different scale 
ranges on angle to the local mean magnetic field direction. 

(|Boldvrevll200f3 ). We will show in the next section that the 
data, in general, support a value closer to —5/3. 

The two decade wide range of scales between kpi — 
3 x 1CT 3 and kpi = 1 X 1CT 1 over which the spectral indices 
vary only slightly can be considered by definition to be the 
inertial range. At the largest values of kpi (blue circles in 
Figure [2]), the spectral index decreases compared to the in- 
ertial range for all 0b > 10° and a value close to —2 is found 
for a wider range of 0b- This is the beginning of the ef- 
fect of kinetic plasma phy sics close to the ion gyroscale (e.g. 
ISchekochihin et al. I [2003 ') . as we show in the next section. 
Thus, we have covered the entire inertial range of MHD tur- 
bulence, from the effective outer scale defined by the break 
point from 1/f scaling at the lowest frequencies, to the be- 
ginning of kinetic effects at the ion gyroscale. 



4 DEPENDENCE ON HELIOCENTRIC 
DISTANCE 

As the solar wind travels away from the Sun, the ion gy- 
roscale increases in size, the mean magnetic field strength 
decreases, the collisional and turbulent age of the plasma 
increases, and the outer scale moves to lower frequencies. 
Here we perform the same analysis as in the previous sec- 
tion for five 50 day periods of Ulysses data from 1995, taken 
at different distances from the Sun and so investigate the 
dependence on heliocentric distance of the anisotropy and 
spectral break points. 

A summary of mean plasma parameters for each pe- 
riod is given in Table \T\ as well as the anisotropic spec- 
tral indices and the mean of the normalised cross helicity, 
a c = (2v- b)/(v 2 +b 2 ) where v and b are fluctuations of the 
velocity and magnetic field in Alfven units; a c is calculated 
spectrally for the whole of each period and the values quoted 
here are the mean values between 6 x 10~ 5 and 8 x 10 Hz, 
which corresponds to the lowest frequency wavelet used. The 
values of a c are consistently above zero, showing the im- 
balanced nature of the turbulence. The values of the spec- 



Figure 3. Dependence of power anisotropy for each of the five 
periods in Table [T] on wavenumber scaled to the ion gyroscale, 
kpi. 



tral indices in Table [T] calculated in a fixed range of scales 
2 x 10~ 2 ^ kpt ^ 2 x 10 _1 for all five periods, are similar and 
broadly consistent with —5/3 and —2 found in the previous 
section. Due to the changing gyroradius and magnitude of 
B, this would not be true if we used a fixed range of frequen- 
cies, we would instead see increasingly negative values of the 
spectral index as the solar wind moves away from the Sun. 
This is because as the ion gyroradius increases with distance 
from the Sun a set frequency can become comparable to, or 
smaller than the ion gyroscale, and thus slip i nto the dis- 
sipatio n range, where the spectra are steeper (|Chen et alj 
l2010bh . 

Figure [3] shows the power an i sotrop y, Pi / P\\, a quan- 
tity first introduced by IPodestal (|2009l ). as a function of 
kpi for the five periods. The data from different periods 
collapse onto a single line, particularly at larger wavenum- 
bers. Its slope at scales within the inertial range is approx- 
imately 1/3, which is the value expected for a critically 
balanced turbulence of lGoldreich fc Sridharl l|l995h . because 
k~ 5/3 /k~ 2 = fc 1/3 . There is a peak at kp i = 0.4 and a tr ough 
at kpi = 0.7, consistent with Figure 6 of lPodestal {uw.j. . For 
kpi < 10 -2 the values of Pj_/P|| mostly rise above the k 1 ^ 3 
line and fill the region 1 < Pj_/P|| < 2 which is the region 
of near-isotropy corresponding to the roll over to a spectral 
index of — 1 shown in Figure [1] for the closest period to the 
Sun. 

All five periods have a similar peak power anisotropy 
of Pl/P m ~ 5 . This is close to the value of 7 reported by 
IPodestal (|2009D for power aniso tropy in t he fa st wind in the 
ecliptic at 1 AU. We note that IPodestal (|2009l ) used smaller 
angular bins than we have here, which could decrease the 
measured power in the P|| angle bin, thus increasing the 
anisotropy (there is likely to be little difference in P± due to 
the approximate invariance of power at angles close to 90°). 

We see that the peak power anisotropies of all five pe- 
riods are very similar. If we assume a universal anisotropic 
scaling in the inertial range and take the result shown above 
that the outer scale is nearly isotropic, these results im- 
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Start R \B\ Vsw Va Ui &c Ti Pi Inertial Range a 

Day (AU) (nT) (km s" 1 ) (km s" 1 ) (cm" 3 ) (xl0 5 K) (xl0 3 km) a(P ± ) a(Pj|) 



100 


1.48 


2.82 


758 


56 


1.21 


0.65 


2.49 


1.52 


-1.62 ±0.03 


-1.92 ±0.02 


150 


1.76 


1.98 


780 


19 


0.77 


0.66 


2.09 


1.98 


-1.58 ±0.02 


-2.00 ±0.01 


200 


2.11 


1.46 


785 


45 


0.50 


0.57 


1.91 


2.57 


-1.64 ±0.01 


-1.9 ±0.1 


250 


2.46 


1.17 


779 


41 


0.38 


0.52 


1.75 


3.06 


-1.69 ±0.01 


-1.9 ±0.1 


300 


2.77 


0.97 


766 


38 


0.31 


0.46 


1.61 


3.56 


-1.64 ±0.02 


-1.94 ±0.04 



Table 1. Results from five different periods of Ulysses data from 1995, showing start day and mean physical parameters for each period. 
The spectral indices of the parallel and perpendicular inertial range arc calculated fitting straight lines to the logarithm of the spectra 
in the range 2 X 10~ 2 < kpi < 2 X 10 _1 . 



ply, perhaps surprisingly, that the total width of the inertial 
range does not change significantly with distance from the 
Sun. Indeed Figure [3] shows that the range of scales from 
where Pj_/Pu = 1 to kpt = 1 is always approximately two 
decades (5 x 1(T 3 < k Pi < 5 x 10" 1 ) Wh i le (as previous 
studi es have shown (|Bavassano et all 1 19821 ; iHorburv et al.l 
Il996f ) both the outer scale and the ion gyroscale increase 
with distance from the Sun (Table [T]), they increase in such 
a way as to keep the total width of the inertial range ap- 
proximately constant. 

Let us make a simple estimate of how the ratio of pi 
and the outer scale (L) varies with distance from the Sun 
R. We assume |S| oc i? -148 , |V| w constant, T oc R~ 102 
and L oc R 1 , the sca lings that have be e n obtained from 
Ulysses observations |Ebert et al.l 12009] ; iGoldstein et al.l 
ll99fj ; lHorburv et al.lll99fj ). Then: 

L oc W. oc (2) 

This very weak dependence on R is essentially unmeasurable 
due to the scatter in the power anisotropy measurements at 
large scales and the small range of heliocentric distances cov- 
ered. This explains why we do not see a significant increase 
in the width of the inertial range. 

This simple analysis taken in conjunction with the scal- 
ing of the magnetic field strength from the Parker spiral 
equation also leads to the conjecture that the inertial range 
might be expected to be wider closer to the Sun. Close to 
the Sun the magnetic field decreases like \B\ oc R~ 2 , whereas 
further out in the h eliosphere it decreases with \B\ oc i? -1 
jBurlagall 19841 , 120021 ). This would imply L/pt oc R~ 1/2 close 
to the Sun and L/ pi oc R 1 ^ 2 further out in the heliosphere. 
The Ulysses results appear to be in the transition region be- 
tween these two behaviours where the L/ P i ~ const. Thus, 
the inertial range in the corona could be wider than that 
observed in fast solar wind at 1 AU and wider again in the 
outer heliosphere. The winding of the Parker spiral controls 
the scaling of the magnetic field magnitude with radius and 
so this scaling is also dependent on heliospheric latitude. 

Finally, we show the similarity between the five periods 
we analyse and the effect of power enhancement parallel to 
the field close to the ion gyroradius by plotting compensated 
spectra for the parallel and perpendicular power. Compen- 
sated spectra are defined: 



Pc(k) = 



P{k) ( k_ 
ko 



(3) 



where kopi = 2x 10 and P±o = P±(koPi); indicated in Fig- 
ure [4] by an arrow. The spectral indices used to compensate 
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Figure 4. Perpendicular and parallel power for each of the five 
periods in Table [l] compensated to remove a spectral gradient of 
—5/3 from the perpendicular power and —2 from the parallel. 



the spectra are a — —5/3 for the perpendicular spectrum 
and a — — 2 for the parallel spectrum. 



Both P± and Pu show a horizontal region between 



2x 10" 



< 



kpi < 5 x 10 1 implying that the spectral indices 
-2 for the perpendicular and paral- 



are close to —5/3 and 
lei power, respectively. At kpi < 10~ 2 the parallel spec- 
trum turns downward, the steeper gradient indicating the 
roll over towards isotropy and a spectral index of —1 has 
begun. The outer scale break point is at kpi w 10~' ! where 
another downward turn is seen in both the perpendicular 
and parallel spectra. 

At kpi ~ 0.7, a clear peak can be seen in the parallel 
power which coincides with the steepening of the perpendic- 
ular spectrum. This is remarkably consistent across all five 
periods used. The parallel power peak suggests a local en- 
hancement (perhaps injection) of energy possibly due to ion 
kinetic instabilities, but does not of course prove it. Ulysses 
high speed wind particle distr ibutions are often clo se to the 
firehose instabi lity threshold jMatteini et al.ll2007l . see also 
simulations by iHellinger fc Travnicekll2008l ). The enhance- 
ment in P|| at kpi m 0.7 is consistent with fluctuations 
whose wavevectors are parallel to the field, although the 
larger power when 6b ~ 90° means that such an enhance- 
ment would not be detectable in the perpendicular power, 
so the total wavevectors may be oblique. What is clear is 
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that the anisotropy is reduced at this scale, perhaps by an 
admixture of fluctuations with kit ^ k±. 

Note that the peak in the parallel power is clearly in the 
same place for all spectra plotted using the gyroscale nor- 
malisation, despite the fact that frequencies and power levels 
change significantly over the 5 periods used. This means that 
what we see is a physical effect associated with the gyroscale 
and not an instrumental or spacecraft effect. 



5 CONCLUSIONS 

We have used a wavelet analysis of Ulysses fast polar so- 
lar wind observations to study the anisotropy of the entire 
inertial range of turbulence with respect to the local mean 
magnetic field. At the outer scale turbulence is isotropic. The 
spectral index is approximately —1 for scales larger than the 
outer scale and changes to —5/3 at smaller scales. 

At scales smaller than the outer scale, the turbulence 
develops a spatial anisotropy with smaller power in fluc- 
tuations that vary along the local mean magnetic field 
and a corresponding spectral index of — 2, while the larger 
power in cross-field fluctuations has a —5/3 spectrum. These 
spectral scalings are consistent with a criti cally balanced 
Alfvenic cascade (|Goldreich fc Sridharlll995h . The perpen- 
dicular scaling appears close r to —5/3 than to —3/2 advo- 
cated by more re cent theory (|Boldvrevll2006l ) and numerical 
simulations (e.g. iMaron fc Goldreichll200lf ). 

The width of the inertial range remains constant at all 
heliocentric distances investigated. This can be understood 
by considering the change in the ratio of outer scale to the 
ion gyroscale with heliocentric distance and leads to the in- 
teresting conjecture that the inertial range is wider closer to 
the Sun and in the outer heliosphere (see Section 4). 

The small scale end of the inertial range seems natu- 
rally to scale with ion gyroradius. There appears to be an 
injection of power into parallel wavenumbers at kpi ~ 1, 
conceivably due to the firehose instability. 

Our results show yet again that the solar wind is a 
unique laboratory for studying plasma turbulence across a 
wide range of scales in extraordinary detail. The isotropy 
at the outer scale and the remarkable consistency of the 
behaviour at the ion gyroscale show that solar wind turbu- 
lence, perhaps surprisingly, is quite a 'clean' case and so is 
not only interesting for its own sake, but also as a represen- 
tative case study of a magnetised plasma of the kind that is 
ubiquitous in the Universe. 

We would like to thank Miriam Forman, Petr Hellinger 
and Lorenzo Matteini for useful discussions. This work was 
funded by the STFC and by the Leverhulme Trust Network 
for Magnetized Plasma Turbulence. 



REFERENCES 

Alexandrova O., Janssen J., Carbone V., Veltri P., Sorriso- 

Valvo L., 2008, ApJ, 674, 1153 
Balogh A., Beek T. J., Forsyth R. J., Hedgecock P. C, Mar- 

quedant R. J., Smith E. J., Southwood D. J., Tsurutani 

B. T., 1992, A&AS, 92, 221 
Bavassano B., Dobrowolny M., Mariani F., Ness N. F., 

1982, J. Geophys. Res., 87, A5, 3617 



Bieber J., Wanner W., Matthaeus W. H., 1996, J. Geophys. 

Res., 101, 2511 
Boldyrev S., 2006, Phys. Rev. Lett., 96, 115002. 
Burlaga L. F., Klein L. W., Lepping R. P., Behannon K. 

W., 1984, J. Geophys. Res., 89, A12, 10,659 
Burlaga L. F., Ness N. F., Wang Y.-M., Sheeley Jr. N. R., 

2002, J. Geophys. Res., 107, All, 1410. 
Chapman S. C, Hnat B., 2007, Geophys. Res. Lett., 34, 

L17103 

Chen C. H. K., Wicks R. T., HorburyT. S. , Schekochihin 

A. A., 2010, ApJ, in press, preprint: larXiv:0909.2683V l 
Chen C. H. K., Horbury T. S., Schekochihin A. A., Wicks 

R. T., Alexandrova O., Mitchell J., 2010, (submitted), 

Phys. Rev. Lett., ? 
Cho J., Lazarian A., 2004, ApJ, 615, L41 
Cho J., Lazarian A., 2009, ApJ, 701, 236 
Dasso S., Milano L. J., Matthaeus W. H., Smith C. W., 

2005, ApJ, 635, L181 
Ebert R. W., McComas D. J., Elliott H. A., Forsyth R. J., 

Gosling J. T., 2009, J. Geophys. Res., 114, A01109 
Fredricks R. W., Coroniti F. V., 1976, J. Geophys. Res., 

81, 31, 5591 
Goldreich P., Sridhar S., 1995, ApJ, 438, 763 
Goldreich P., Sridhar S., 1997, ApJ, 485, 680 
Goldstein M. L., Roberts D. A., Matthaeus W. H., 1995, 

ARA&A, 33, 283 
Goldstein B. E., Neugebauer M., Phillips J. L., Bame S., 

Gosling J. T., McComas D., Wang Y.-M., Sheeley N. R., 

Suess S. T., 1996, A&A, 316, 296 
Hellinger P., Travmcek P. M., 2008, J. Geophys. Res., 113, 

A10, A10109 

Horbury T. S., Balogh A., Forsyth R. J., Smith E. J., 1996, 

A&A, 316, 333 
Horbury T. S., Forman M. A., Oughton S., 2008, Phys. 

Rev. Lett., 101, 175005 
Lithwick Y., Goldreich P., Sridhar S., 2007, ApJ, 655, 269 
Maron J., Goldreich P., 2001, ApJ, 554, 1175 
Matteini L., Landi S., Hellinger P., Pantellini F., Maksi- 

movic M., Velli M., Goldstein B. E., Marsch E., 2007, 

Geophys. Res. Lett., 34, L20105 
Matthaeus W. H, Oughton S., Pontius Jr. D. H., Zhou Y., 

1994, J. Geophys. Res., 99, AI0, 19,267 
Milano L. J., Dasso S., Matthaeus W. H., Smith C.W., 

2004, Phys. Rev. Lett., 93(15), 155005 
Nicol R. M., Chapman S. O, Dendy R. O., 2009, ApJ, 703, 

2138 

Osman K. T., Horbury T. S., 2007, ApJ, 654, L103 
Podesta J. J., 2009, ApJ., 698, 986 

Sahraoui F., Goldstein M. L., Robert P., Khotyaintsev Y. 

V., 2009, Phys. Rev. Lett., 102, 231102 
Schekochihin A. A., Cowley S. O, Dorland W., 2007, 

Plasma Phys. Control. Fusion, 49, A195 
Schekochihin A. A., Cowley S. O, Dorland W., Hammett 

G. W., Howes G. G., Quataert E., Tatsuno T., 2009, 

ApJS, 182, 310 
Taylor G. I., 1938, Proc. Roy. Soc, 164, 919, 476 
Torrence O, Compo G. P., 1998, Bulletin of the American 

Meteorological Society, 79, 1, 61 
Weygand J. M., Matthaeus W. H., Dasso S., Kivelson M. 

G., Kistler L. M., Mouikis O, 2009, J. Geophys. Res., 114, 

A07213 



